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Abstract. Carey’s Equality pertaining to stationary models is well known. In this paper, we have 
stated and proved a fundamental theorem related to the formation of this Equality. This theorem 
will provide an in-depth understanding of the role of each captive subject, and their corresponding 
follow-up duration in a stationary population. We have demonstrated a numerical example of a 
captive cohort and the survival pattern of medfly populations. These results can be adopted to 
understand age-structure and aging process in stationary and non-stationary population population 
models. Key words: Captive cohort, life expectancy, symmetric patterns. 


1. Introduction 

The motivation to explore the question from which Carey’s Equality ultimately emerged stemmed 
from the lack of methods for estimating age structure in insect populations in particular, and in a 
wide range of animal populations generally. Although various methods exist for estimating the age 
of individual insects, including the use of mechanical damage, chemical analysis and gene expres¬ 
sion, many are costly, most require major training and calibration efforts and none are accurate 
at older ages. Because of these technical constraints, our concept was to explore the possibility 
of using the distribution of remaining lifespans of captured insects (e.g. fruit flies) to estimate 
age structure statics and dynamics. The underlying idea was that the death distribution of all 
standing populations must necessarily be related to the population from which it is derived, ceteris 
perabus. 

James Carey first observed the symmetric survival patterns of a captive cohort and follow-up 
cohort of medflies (see for example in |21 H|). Later it was referred to as Carey’s Equality [3] 
although a formal mathematical statement on a generalized Equality of Carey’s type was missing. 
The phenomena of Equality of mean length of life lived by a cohort of individuals up to a certain 
time, and the mean of their remaining length of life in a stationary population, was well established 
mm- Proof of Carey’s Equality in [3] is based on the stationary population property explained in 
m although similar phenomena were also found to be useful in understanding the renewal process, 
see Chapter 5 in [13]. We have formally stated a fundamental theorem to justify Carey’s Equality 
under a stationary population model, which would allow us to obtain expected lifetimes of a captive 
cohort of subjects in which each subject is captured at a different age, starting at birth. Further, we 
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have proved our statement under a certain combinatorial and stationary population modeling set¬ 
up. These results are extended to two dimensions for estimating age-structure of wild populations 
and understanding internal structure of aging process, age-structure of wild populations [3]. One 
also needs to understand the implications of our results for estimating the age structure in stable 
populations [5] . 


2. Main Theorem 

Theorem 1. Suppose ( X , Y. Z ) is a triplet of column vectors, where X = [xi, X 2 , ■ ■ ■ , Xk] T , Y = 
[?/i, 2 / 25 ' * • ,yk\ T , Z = [zi, Z 2 , ■ • • ,Zk] T representing capture ages, follow-up durations, and lengths 
of lives for k—subjects, respectively. Suppose, F(Z), the distribution function of Z is known 
and follows a stationary population. Let G\ be the graph connecting the co-ordinates of Sy, the 
survival function whose domain is N(k) = {1, 2, 3, k} i.e. the set of first k positive integers and 
Sy (j) = Vj for j = 1,2,..., k. Let G 2 be the graph connecting the co-ordinates of Cx , the function of 
capture ages whose domain is N(k ) and Cx ( j ) = Xj for all j = 1,2,..., k. Suppose C x (— j ) = Xj 
for all j = 1,2 ,...,k. Let Li be the family of graphs constructed using the co-ordinates of C* x 
consisting of each of the k\ permutations of graphs. Then one of the members of Li ( say,H g ) is a 
vertical mirror image of G\. 

Proof. In the hypotheses, the information on capture ages, follow-up durations, and lengths of 
lives are given as column vectors. For analyzing data collected from life table studies on different 
organisms (including life expectancy), column vectors are used for representing a convenient data 
structure, where the matrix of interest has only one column. Capture ages of all subjects in our 
paper are arranged in a column such that it will be suitable to do addition operations on these 
ages with a column vector consisting of lives left for corresponding subjects. One also can compute 
inner product of these two column vectors, but the resultant scalar obtained by such operation is 
not of interest in the present context. Life table analysis in demography considers column-wise 
information and each column has a special purpose even though the information obtained from the 
last column has the key data needed for actuarial modeling, mortality analysis, etc, Using these 
column vectors of capture ages and follow-up durations from the hypothesis, we have constructed 
two sets, 1 1 and J 2 . 


GENERALIZATION OF CAREY’S EQUALITY AND A THEOREM ON STATIONARY POPULATION 4 
Let Ii = {(i,yi) : 1 < i < k} and J 2 = {(z, #j) : 1 < i < k} . G\ is constructed using specific 
ordered pairs of ii, explained later in the proof, and (SL is constructed by corresponding ordered 
pairs of J 2 . There are two criteria, U\ and f/ 2 , that govern the one-to-one correspondence properties 
of the two functions, S and Cx- 


U\ = {No two subjects are captured who are of the same age and 

there are no two subjects whose follow-up times until death are identical} . 


Ui = {There is more than one subject which has the same capture ages, and 

there is more than one subject which has the same follow-up durations until death} . 

Suppose fi and / 2 are probability density functions of capture times and follow-up times to 
death, then by [3], we will have fi(a) = / 2 (a), which is the probability of an individual who lived 
'a' years is the same as probability that this individual lives 'a' years during the follow-up |3j. For 
an infinite population in a continuous time set-up for each x G Cx there exists a y E Sy such that 


( 2 . 1 ) 


X 


[s)ds + 




z(s)ds 


and 


1 r°° 

(2.2) x A =y A = - J z(s)ds, 


where x A and y A are the average capture ages and average follow-up durations until death. Sum 
of the age at capture of i th — subject (i.e. Xj) and the follow-up duration (or time left in the current 
context) of i th —subject (i.e. yi) is equal to the total length of the life (i.e. Zi) for the i th — subject. 
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See [31 [6] and Chapter 3 in [?j for a description of stationary population models and see page 49 
in [Sj and see page 52 in [9] for the stationary distributions. 

Let’s define max = yt x for some t\ h — subject out of k— subjects, 

max {y t : for all i G N(k) and i ^ H} = y t2 for some t^—i subject out of k — 1 subjects and so on 
until we arrive at, 

max [y t : for all i G N(k ) and i ± h,t 2 , t(fc-i)} = Vt k - 

G i is constructed by joining the co-ordinates of the set S on the first quadrant, where 


(2.3) 


S = { (ti.Vh), ft;, 


We call each co-ordinate of S' as a cell of S. One can visualize, cells in S are made up of ordered 
pair of co-ordinates, where abscissa is the subject captured and ordinate is the life left for this 
subject after capture. Here, the graph, we mean by a curve obtained by joining the cells in S. 
Each cell, except the first and the last cell of S are joined to both sides of its neighboring cells. 
In case there is more than one subject with a maximum value at one or more stage above, it will 
lead to two or more identical co-ordinates that are used in G\. When S satisfies Ui, G\ is a graph 
of decreasing function; when S is satisfies criterion f/ 2 , G\ is a graph of combination of decreasing 
and non-decreasing functions. 

We can construct a sequence of quantities x ui ,x U2 , • • • ,x Uk similarly to S to form the set T, given 
below: 


(2.4) T = {{ui,x Ul ) , (u 2 ,x U2 ) , ■ ■ ■ , (u k ,x Uk )} 

Corresponding to (H, y tl ) G S there exists a (j, x 3 ) G I 2 which could be (u\, x Ul ) or some other cell 
in T. Note that, Cx = {(— j, Xj) : 1 < j < &}. Corresponding to (ti,y tl ) G S if there exists a cell 
(|— i\ ,Xi) G such that |H — |—*|| = 0 and | y tl — Xj| = 0 then we call and (|— i\ ,Xi) are a 

pair of equidistant cells from a vertical axis. We can construct a graph Hi using cells of C* x in the 
natural order of integers. In total, we can have k\ permutations of orders to construct Hi ... Hk\- 
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Suppose we denote the first combination of cells as j y— j^\ xfj : 1 < j < kj and the second 
combination of cells as j j( 2 \ x^^j : 1 < j < A; j and so on. Negative index here indicates that 
the subject captured is considered on the negative x— axis, which is similar to the left part of the 
Figure 15. 11 if we consider the values 25, 50 75 and 100 for pre-capture segment as -25, -50, -75 
and -100 as we visualize them on the negative x— axis. Thus by previous arguments, a family of 
graphs "H is constructed. One of these combinations, for example, the g th combination, is used to 
construct a graph which we denote with H g: which satisfies following Equality: 


ti- |-1 

t'2 — I —2 


( 9 ) 

(9) 


= 0 and 
= 0 and 


(. g ) 

yt! - x\ y> 

( g ) 

yt 2 x 2 


= 0 
= 0 


4-1-1 


( 9 ) 


= 0 and 


(9) 

Vt k ~ x k 


= 0 


here j ^|— 1|^ , , ^|— 2\^ , — k\^ , | are g th combination of cells such that 

H g is a vertical mirror image of G\. Image of G\ is visualized as if it is seen from the mirror kept 
on y- axis. Note that, we have generated k\ graphs by our construction, and one of such graph, 
which we called as H g is shown to have vertical mirror image of G\. □ 


Our theorem establishes existence of graph depicting mirror images of the pre- and post-capture 
longevity distributions in Figure l57ll In general, in the area of population biology where survivor¬ 
ship curves of captive cohort obey Carey’s Equality, there is a need for understanding the pattern 
within the data structure. The results by 13 HU explain functional symmetries as an application of 
renewal theory. We do not depend on any of the classical works on renewal equation and renewal 
theory proposed by Lotka TB, [Da] [T7j and by Feller [18] (also see Chapter XI in [IS]). Our in¬ 
spiration is purely from experimental observations demonstrated by James Carey and a statement 
on stationary populations in the equation (1) in the paper by Vaupel [3], However, our theory 
and method of proof uses sequentially arranged data of captive individuals, which was also usually 
done in renewal theory analysis or proving renewal-type of equations. 
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Hence our method provide an alternative and independent approach for such kind of sequentially 
arranged captive subjects. Besides relating the captive ages and corresponding follow-up durations 
of subjects in a stationary population, the principles of the theorem helped to visualize person- 
years in a follow-up starting at birth, in a stationary population model in which subjects of each 
age are captured. 


2.1. Life expectancy. First we construct the age structured survival function using l(s), the 
number of captured subjects at age s and Z(so), the number of subjects at the beginning of the 
follow-up. Suppose So is the time at the beginning of the follow-up, s, is the time at the i th time 
point of observation for i — 1 , 2 , • • • , k. Suppose each of the y t . for i = 1 , 2 , • • • , k fall exactly in 
one of the time intervals (sj_i,Sj), then the expectancy of life is pp. The probability of death 
over the time period are g(s*) = pp with the survival pattern, Z(s;) = k — i for i — 1 , 2 , • ■ ■ , k — 1 
and Z(sfc) = 0, q(sk ) = 1. Suppose there are n* number of y ti falling within (sj_i,Sj) such that 
'Y^=i n i = k and Z(sj) follows the previous construction. If n.j > 1 in one or more of the (sj_i,Sj) 
then there must exist empty cells where the event of death is avoided. Let us define a number c\ 
as follows: 


Ci = {Number of cells (sj_i, Sj) for i — 1,2, • • • , k where exactly one of the y ti falls } 


If Ci = k, then the life expectancy is pp. If ci 7 ^ k, then at least one of the (sj_i,Sj) is empty. 
There could be several combinations of distribution of y tl in (sj_i, Sj) when the event Ci 7 ^ k occurs, 
and we explain in the following remark one possible situation in which deaths are concentrated in 
the early ages and at late ages. Other combinations can be evaluated using similar constructions. 


Remark 1. Suppose n\ — n 2 = • • • = n ri — 1; n ri+ 1 = n ri+ 2 = • • • = n r2 = 0; and n r2+ 1 = n r2+ 2 = 
• • • = rik = 1 such that X^=i n * = ^ f° r some 1 < T]_ < r 2 < k. The probability of death at various 
ages, q(si), is, 
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k—i 


for i — 0,1, ■ ■ • , (rr — 1) 


(2.5) 


q(si) 


0 fori = n, (rr + 1), • • • , (r 2 — 1) 

< 


i 

k—(ri+i) 


for i — 0,1, • • • (k — 2) 


1 for i = (k — 1) 


The life expectancy, e(so), is obtained by the formula below: 


e(s 0 ) = 


2k 


2ri — 1 2fc—1 

i 2 k ~i}+ { 2k - ft 

j=2ri+l 


i=l 


+ 


(k - ri) (r 2 — 7T — 1) 
k 


( 2 . 6 ) 


2.2. Person-years and means. Under the above set-up, in this section we derive the mean age 
of the captive cohort in terms of the mean of the person-units followed. Suppose a(k. x ) denotes 
k tll —subject (k > 0) is captured at age x, then a(k , x)dk is number of subjects captured at age 
x. The mean age at capture for all subjects of the cohort formed of subjects of all ages, c(0), is 


(2.7) 


c(0) = 


/o°° [ x /i°° a (k, x)dk\ dx 
/o°° fi° a (^’ x )dkdx 


Let /q b(y)dy be the number of deaths out of a(k,x)dkdx during the age 0 to s. The 

number of subjects surviving at age s is 



a(k,x)dkdx 



b(y)dy 
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and the number of subjects surviving at age (n + l)s for some positive integer n is 


r»oo poo 


’0 J 1 


00 /»(n+l)s 

i(k, x)dkdx — £ b(y)dy. 

n=0 -> ns 


Thus, life expectancy at the formation stage of a captive cohort, say, i?[c(0)] can be computed 
by the formula 


( 2 . 8 ) 


£[c(0)] = 


r r - r /, 

/o°° fi° a (^’ x )dkdx 


-(n+l)s 


ds 


We state a theorem for the total person-years to be lived by all the subjects in a birth cohort. 


Theorem 2. Suppose subjects of each age of life in a population are captured. Then, using the 
constructions in c(0) and h?[c(0)] 7 the total person-years, say, T (a, c(0), i?[c(0)]) ; that will be lived 
by newly born subjects in a stationary population can be expressed as 

/ oo /»oo 

J a(k, x)dkdx. 


Proof. We have, 
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"OO />oo 


(c(0) + i±7[c(0)]) / / a(k,x)dkdx = 

Jo J 1 Jo 


x / a(k,x)dk 


dx + 


"OO /»oo /»oo 


'0 JO J1 


00 /" °° r(n+l)s 

a(/c, x)dkdxds — I IEI b(y)dy | ote 

n =0 ' 


'0 \„—nJns 


x J a(k,x)dk 


dx + 


(2.9) 


"OO /*oo 


r 0 Jl 


a(fc, x)dkdx — 


°° /»(n+l)s 


b(y)dy 


n =0 ' 


<is 


The R.H.S. of (12.91) is sum of ages of subjects of all ages in a captive cohort, and total person- 
years to be lived by the captive cohort, which is T (a, c(0), i?[c(0)]). □ 


3. History and Related Results 

An Equality arising out of symmetries of life lived and life left was named as Carey’s Equality by 
James Vaupel f3], to highlight the discovery of certain symmetries in his biodemographic experi¬ 
ments by James Carey (see, for example, mm)- Vaupel [3] and Goldstein HH have proved equality 
on such symmetries as a direct application of renewal theory. Our main theorem in this article 
is not inspired by renewal theory, but we conceptualized our approach directly from experimen¬ 
tal results shown by James Carey and then used equation (1) from Vaupel [3j in our hypothesis. 
Renewal theory has long history even before the seminal works on population dynamics by Alfred 
Lotka (see, for example, [14] [15, H]), who has used an integral equation of type (13.11) to link 
number of births at time t > 0 with number of births a women has at time t = 0. 

(3.1) B(t)= [ G(t,a)da + H(t) 

Jo 

where B(t) is total number of births at time t, G(t,a ) is number of births from a women who 
is at age a and alive at time t and H(t) is number of births from a women who is alive at 1 = 0. 
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Usually, we compute B(t) from a = a to a = /?, where a is lower reproductive age and /? is upper 
reproductive age. G(t, a) can be written as, 


(3.2) G(t, a) = B(t — a)s(a) 

where B(t — a) is total number of births at time (t — a ) and s(a) is chance of surviving to exact 
age a. 13.11 is also referred as renewal equation. Combining 13.11 and 13.21 we can write B(t) within 
reproductive ages as, 


rP 

(3.3) B(t) = / B(t — a)s(a)da 

J a 

William Feller [T9] provided foundations of renewal theory in his book (see Chapters VI and XI) 
as did authors of other books which were written exclusively on renewal theory (see for example, 
|13j ) or contained chapters devoted to basic renewal theory (for example, see Chapter 6 in [9] and 
Chapter 12 in |20|). Several applications of convolutions of independent random variables which 
we see in renewal theory can also be found in understanding disease progression between one stage 
to another stage, epidemic prediction and so forth (for example, see [2T] and [22]). 

4. Example and visualization 

We provide here a practical application based on the medfly population with a visual depiction 
of Carey’s Equality (Figure f5TTT) . Note the symmetry of the pre- and post-capture segments of the 
lifespans of individuals in the population that underlies the equivalency of life-days which, in turn, 
underlies the equivalency noted by |3], "If an individual is chosen at random from a stationary 
population...then the probability the individual is one who has lived a years equals the probability 
the individuals is one who has that number of years left to live.” 

Carey’s Equality is important because it both builds on and complements the properties of 
one of the most important models in demography: the stationary population model. Station¬ 
ary population model is fundamental to formal demography because it is a special case of the 
stable population model, and provides explicit expressions that connect the major demographic 
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parameters to one another, including life expectancy, birth rates, death rates and age structure, 
see Chapter 3 in [7]. A graphical depiction of Carey’s Equality, shown in Figure [5721 shows the 
interconnectedness of these parameters. Note that the shape of the stationary population 15.2( a) 
as well as both its age and death distributions 15.2( b) are identical, and that the proportion of 
each age class in the captive population 15.2( c) is identical to the proportions within the whole 
population (i.e. due to the symmetrical distributions of pre- and post-capture lifespans in 15.2( b)). 

5. Discussion 

The approach we used in this paper to demonstrate the mathematical identity underlying Carey’s 
Equality is fundamentally different than that used by previous authors. Originally James Carey 
created a simple life table model (Table 1 in |2J) to demonstrate the equality of age structure 
and post-capture deaths in a stationary population, the identity of which was then formalized by 
statisticians Hans Muller and Jane-Ling Wang [2], Vaupel |3] and Goldstein mi followed by using 
mathematical first principles to derive the equality. 

Our approach differs from the one we just outlined. Instead of first formulating and then deriving 
the equality, we justified the constancy of post-capture patterns of death in stationary medfly 
populations in the following two steps. The first was to conjecture that in stationary populations 
the ordered pre- and post-capture life course segments of individuals will be symmetrical (Theorem 
[T]) and that total person-years in a captive cohort can be formulated (Theorem [2]). The second 
was to prove these relationships using a series of mathematical assumptions (i.e. the respective 
proofs). 

Our approach contributes to the demographic literature in general and specifically to an under¬ 
standing of Carey’s Equality in several ways. First, our theorems provide an independent method 
for formulating a mathematical relationship in formal demography. Indeed we are unaware of any 
other models in formal demography that have involved proofs from mathematical conjectures (the¬ 
orems). Second, our proofs allow us to state unequivocally that the Carey Equality will be true in 
all stationary populations. Although this is a logical outcome from all of the earlier approaches, 
our approach makes this result both explicit and conclusive. Third, our proofs required that we 
draw on set theory, an area of mathematics involving logic that is not commonly used in demog¬ 
raphy. As a consequence of the problem framed in fundamental mathematics, Carey’s Equality is 
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Days 

FIGURE 5.1. Schematics of a randomly sampled hypothetical stationary Mediter¬ 
ranean fruit fly population (adapted from [ID]). Horizontal lines depict the life 
course of individual flies from birth (eclosion) to death, (a) Example of the life 
courses of captive individuals divided into pre- and post-capture segments and or¬ 
dered from shortest-to-longest (top-to-bottom) post-capture lifespans. Inset shows 
the stationary medfly population in which the yellow tick marks depict the ages at 
which individual medflies are sampled, (b) Same as (a) except with the pre- and 
post-capture segments of individuals decoupled and both ordered according to length 
from top-to-bottom to show the symmetry (mirror image) of the distributions. 

better situated both to draw from and be extended into other areas of basic mathematical theory. 
Fourth, using the ideas of the main theorem, we are positioned to obtain further results related 
to Carey’s Equality such as for higher dimensions and probabilistic and deterministic results for 
multiple captive cohorts. 


6. Conclusions 

Our paper offers new sets of tools, techniques and theoretical framework in terms of visualization 
of the demographic data involving capture ages and development of new theoretical ideas for 
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FIGURE 5.2. Structural interconnectedness of an hypothetical stationary Mediter¬ 
ranean fruit fly population showing: (a) days to death by age groups; (b) age and 
death distributions; and (c) captive survival. 

analyzing data obtained by captive cohorts. Such approaches will have applications in other 
demographic situations, for example, understanding aging patterns in a captive cohort data when 
information on lives left is right truncated, projecting various scenarios of demographic transition, 
etc, 

The main result of our paper will be useful in understanding the relationship between average 
lengths of lives of captured, follow-up, and total lengths of the lives in a stationary population. 
Our method of re-structuring the follow-up durations of captive cohort can be adopted also for 
non-stationary populations, which can be used for understanding the internal structures of the 
population with respect to the age at capture. This will enable us to look deeper into the aging 
process of stationary and non-stationary populations. For each captive cohort of subjects there 
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exists an associated, exact configuration of a combination of coordinates of the survival graphs, and 
this association is dynamic. The right combination of coordinates is dependent on the formation 
of the captive cohort. The idea of a proof through formation of symmetric graphs, combined 
with captive age distribution is novel. We have demonstrated the utility of such thinking in 
understanding symmetric patterns formed of a captive cohort and associated follow-up lengths. 
This strategy was also helpful for us in deriving formulae for expectation of life in a stationary 
population, discretely, and also using multiple integrals. The theory explained here can be adopted 
to both human and non-human populations. 
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